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1.  Introduction 

The  goal  of  this  project  was  to  build  on  an  understanding  of  the  connections  bc- 
tween  knot  invariants,  exactly  solvable  statistical  mechanics  models  and  discrete 
dynamical  systems  toward  an  answer  to  the  question  of  how  early  and  robust  ther¬ 
modynamic  behavior  appears  in  lattice  gas  automata.  Preliminary  work  focussed 
on  developing  our  understanding  of  state  models  in  knot  theory  and  their  relation 
to  statistical  mechanics  models;  this  work  is  described  in  [1]. 

Rather  than  considering  the  complicated  and  somewhat  more  difficult  case  of 
lattice  gas  hydrodynamics  we  began  with  a  simpler  model;  a  reversible  cellular  au¬ 
tomaton  in  one  dimension.  This  particular  automaton  has  an  additive  conserved 
quantity  which  Takesue  has  used  to  construct  a  one  dimensional  statistical  me¬ 
chanics  model  [2].  He  observes  in  simulations  that  the  cellular  automaton  displays 
thermodynamic  behavior  which  is  well  described  by  this  one  dimensional  statistical 
mechanics  model.  In  [3]  and  [4]  we  showed  that  there  is  a  naturally  associated  two 
dimensional  statistical  mechanics  model  as  well,  then  examined  the  properties  of 
this  model,  and  suggested  how  these  properties  may  relate  to  the  behavior  of  the 
dyneimical  cellular  automaton.  This  work  is  described  in  the  following  section. 

Having  applied  these  ideas  successfully  to  the  simple  case  of  a  reversible  cellular 
automaton,  we  moved  on  to  the  more  interesting  case  of  lattice  gases.  Since  our 
approach  equates  the  spacetime  evolution  of  a  dynamical  system  with  an  equilibrium 
configuration  of  a  statistical  mechanics  model  in  one  higher  dimension,  a  model  of 
’t  Hooft  for  two  dimensional  spacetime  with  discrete  local  coordinate  invariance 
was  a  natural  inspiration  [5].  In  [6]  we  detail  the  construction  of  a  family  of  one 
dimensional  lattice  gas  automata  based  on  this  model.  This  work  is  summarized  in 
section  3. 
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2.  Cellular  automata 


\^e  focus  on  reversible  models  for  several  reasons:  Microscopically,  one  expects 
physical  theories  to  be  reversible.  Moreover,  with  the  additional  constraint  of  dis¬ 
creteness,  reversibility  implies  that  phase  space  volume  is  conserved  under  the  evo¬ 
lution;  hence  an  analog  of  Liouville’s  Theorem  applies.  Finally,  when  some  iocally 
defined  and  additive  quantity  is  conserved,  the  system  is  similar  to  a  Hamiltonian 
system,  the  context  in  which  questions  of  ergodicity  and  thermodynamic  behavior 
are  usually  studied. 


Motivated  by  these  considerations,  Takesue  has  begun  an  analysis  of  the  class 
of  one  dimensional  reversible  cellular  automata  with  a  nearest  neighbor  rule  [7].  In 
particular,  consider  the  the  one  dimensional  autorriaton  \yith  a  boolean  variable  at 
each  site  evolving  according  to 


where  /  ;  {0, 1}^  — >  {0, 1}  is  given  by 


/(A,  /2,v)  :=  A  +  1/  —  /Ji/  —  2Xu  -f  Xfjv, 


rule  26R  in  the  usual  conventions  [8],  There  is  a  conserved  quantity  defined  by 


Takesue  uses  this  result  in  one  test  of  thermodynamic  behavior  [2]:  .4ny  initial 
configuration  has  some  total  energy  which  determines  an  inverse  temperature  3-  In 
simulations  he  observes  a  canonical  ensemble  at  this  temperature  for  subsystems. 
More  precisely,  the  dynamical  distribution  of  energies  £„  for  a  subsysten^  of  size 
n-flis  found  to  be  close  to  D{En)t~^^''  where  D{E„)  :=  (numhcT  of  configurations 
with  energy  £'„)/4"'''h  This  phenomenon  is  observed  even  for  fairly  small  values  of 
N  and  n.  Similarly,  he  finds  that  by  simulating  a  heat  bath  at  either  end  of  the 
system,  rule  26R  supports  a  temperature  gradient  [9], 

Thus  this  automaton  exhibits  the  type  of  surprising  thermodynamic  behavior 
we  would  like  to  understand:  Despite  being  far  from  the  thermodynamic  limit  (small 
N  and  n)  the  dynamical  system  displays  equilibrium  behavior  well  approximated  by 
this  statistical  mechanics  model.  Moreover,  such  success  suggests  that  the  system 
is  ergodic.  We  know,  however,  that  It  cannot  be  truly  ergodic,  for  there  are  0(2’^ ) 
cycles  [7]  and  only  2-V  energy  surfaces,  which  means  that  the  energy  surfaces  must 
be  partitioned  into  many  orbits. 

Our  suggestion  is  that  the  origin  of  this  thermodynamic  behavior  might  be 
found  in  a  naturally  associated  two  dimensional  statistical  mechanics  model.  Tlie 
Boltzmann  weight  Wf}  :=  exp[-3£(o,  6,  c,  d)]  can  be  identified  as  the  transition 
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matrix  for  the  scattering  process  (crj“^ ,  )  — y  (a^,  )  or,  equally  well,  as  the 

Boltzmann  weight  of  a  vertex  model  on  a  square  lattice.  Thus  the  deterministic 
dynamical  evolution  of  the  cellular  automaton  is  modelled  probabilistically  (hence 
not  exactly)  by  two  such  statistical  mechanics  models  on  a  pair  of  dual  lattices 

Now,  recall  the  transfer  matrix  method  for  the  solution  of  two  dimensional  sta¬ 
tistical  mechanics  models  [10].  For  the  model  we  are  considering  here  the  partition 
function  is 

^  =  E  n 

stales  vertices  i 

where  the  energy  Ei  associated  to  a  vertex  i  is  determined  by  evaluating  the  ex¬ 
pression  for  the  conserved  quantity  above  on  the  spins  (boolean  values)  assigned  to 
the  incident  edges  by  a  given  state.  To  apply  the-transfer  matrix  method  to  this 
model  we  would  proceed  by  defining  the  diagonal-to-diagonal  transfer  matrices  and 
find  that  the  partition  function  again  becomes  the  trace  over  a  product  of  identical 
matrices. 


Thus  the  problem  of  solving  this  model  in  the  thermodynamic  limit  reduces  to 
the  problem  of  determining  the  largest  eigenvalue  of  some  2'"^  x  2-^'  transfer  matrix 
T.  Because  commuting  matrices  have  common  eigenvectors  it  has  proved  sufficient 
in  models  like  this  one  to  recognize  T  as  an  element  in  a  one  parameter  family  T(u:} 
of  mutually  commuting  transfer  matrices:  [T(u;),  r(cj')]  =  0,  where  the  spectral 
parameter  u;  is  determined  by  the  coupling  constants/interaction  strengths  in  the 
energy  functional  defining  the  model.  Further,  it  is  easy  to  see  that  a  sufficient 
condition  for  the  existence  of  commuting  transfer  matrices  is  that  the  Boltzmann 
weights  satisfy  the  Yang-Baixter  (star-triangle)  equation  [10.11]: 


(-■)■ 


The  relevance  of  this  result  to  our  problem  is  the  following.  The  two  com¬ 
ponent  (spin  i)  solutions  to  the  Yang-Baxter  equation  have  been  classified.  The 
most  general  solution  is  the  eight-vertex  model  in  which  the  Boltzmann  weights 
are  parameterized  by  elliptic  theta  functions  [12].  In  this  model  the  temperature 
is  related  to  the  theta  function  nome  p  in  such  a  way  that  p  —  Q  corresponds  to 
T  -  'Tcriticai-  Moreover,  only  at  p  =  0  is  the  limit  a;  — »  oo  defined.  In  this  limit  one 
obtains  the  six-vertex  model  and  the  Yang-Baxter  equation  becomes  independent 
of  the  spectral  parameter: 


This  is  precisely  the  condition  realizing  the  type  III  Reidemeister  move  of  knot 
theory  [11,13]  and  is,  as  we  showed  in  [3],  satisfied  by  the  Boltzmann  weights  of  our 
two  dimensional  statistical  mechanics  model.  Thus  our  cellular  automaton  is  really 
the  (critical)  six-vertex  model  in  disguise. 
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Now,  the  dominant  characteristic  of  criticality  is  scale  or  conformal  invari¬ 
ance,  which  means  that  probability  distributions  of  states  are  invariant  and  correla¬ 
tion  functions  transform  covariantly  under  conformal  transformations.  Thus,  if  our 
model  were  really  at  a  critical  point,  i.e.,  N  — ►  oc  as  well  as  T  —  Tcrhical.  the  local 
distribution  of  states  would  be  that  of  the  thermodynamic  limit,  a  result  consistent 
with  an  observed  canonical  distribution.  For  finite  .V  tViis  will  V>e  moderated  by 
finite-size  effects  which  round  off  the  critical  singularity.  Some  finite-size  effects  can 
be  computed  explicitly  [14];  for  example,  the  correlation  length  will  be  asymptot¬ 
ically  proportional  to  A'.  This,  in  fact,  is  exactly  vhat  Takesue  observes  when  he 
measures  the  autocorrelation  function  of  the  energy  flux  to  check  the  Green-Kubo 
formula  [9].  So  the  criticality  of  the  associated  two  dimensional  statistical  mechan¬ 
ics  model  is  nicely  consistent  with  the  observed  thermodynamic  behavior  of  the 
cellular  automaton. 


3.  Lattice  gases 

Moving  on  to  the  simplest  type  of  lattice  gas  model,  in  [6]  we  described  the  con¬ 
struction  of  a  family  of  one  dimensional  lattice  gas  automata  based  on  a  model  of 
’t  Hooft  [5].  In  an  effort  to  clarify  which  features  of  the  lattice  gas  automata  are  the 
consequences  of  which  constraints,  we  proceeded  systematically:  first  setting  the 
kinematics  of  the  model,  then  exploring  the  effect  of  requiring  that  the  dynamics  be 
local,  and  finally  constructing  the  possible  dynamics  consistent  with  the  imposed 
constraints. 

The  most  interesting  of  these  three  steps  is  the  second.  By  local  dynamics  we 
mean  that  two  particles  interact  only  when  they  coincide  in  both  space  and  time. 
The  simplest  consequence  of  this  requirement  is  illustrated  in  Figure  1.  There 
particles  1  and  2  interact  at  time  t  -f  1,  but  particles  1'  and  3  move  past  each 
other  without  interacting  because  their  spatial  coordinates  always  differ  by  an  odd 
number,  precluding  coincidence.  This  apparently  trivial  observation  is  at  the  origin 
of  the  “spurious”  conserved  quantities  often  found  in  lattice  gas  models;  here  the  gas 
decouples  into  two  noninteracting  gases,  one  occupying  the  even  nodes  of  the  lattice 
and  the  other  the  odd.  In  [3]  and  [4]  ive  used  the  conserved  quantity  of  the  cellular 
automaton  to  define  the  weights  of  the  equivalent  statistical  mechanics  model.  Here, 
rather  than  follow  the  usual  approach  to  dealing  with  a  spurious  conserved  quantity 
(modifying  the  model  in  an  attempt  to  break  the  unwanted  conservation  law)  we 
use  it  to  effect  a  radical  transformation  of  the  model  to  one  which  is  again  equivalent 
to  a  2  dimensional  statistical  mechrmics  model.  The  transformation  is  described  in 
[Cj;  to  our  knowledge  this  is  the  first  time  such  an  approach  has  been  taken. 

The  conclusion  is  that  the  model  lies  on  a  first  order  critical  line  between  dis¬ 
ordered  and  totally  ordered  regimes  of  the  six  vertex  model.  We  observe  that  this 
result  applies  to  ensembles  of  systems,  in  the  thermodynamic  limit.  That  the  cor¬ 
relation  length  vanishes,  for  example,  does  not  mean  that  there  is  no  correlation 
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Figure  1. 


between  states  of  the  probabilistic  lattice  gas  automaton  at  successive  time  steps. 
Nor  would  a  nonvanishing  correlation  length  mean  that  there  were  necessarily  space¬ 
like  correlations  in  a  given  initial  state.  Instead,  the  interpretation  is  that  deep  in 
the  interior  of  a  system,  far  from  the  boundary  conditions  and  thus/cr  subject  to 
essentially  random  boundary  conditions,  a  domain  of  sufficiently  large  size  is  in 
thermodynamic  equilibrium.  Thus  our  results  provide  evidence  for  the  first  step 
in  the  derivation  of  the  macroscopic  equations  of  motion— the  assumption  of  local 
thermodynamic  equilibrium,  as  well  as  demonstrating  the  efficacy  of  constructing 
and  analyzing  lattice  gas  automata  according  to  (spacetime)  symmetry  principles. 

In  ‘t  Hooft’s  original  model  the  particle  trajectories  were  to  be  interpreted 
as  defect  lines  in  a  (spacetime)  lattice  which  developed  curvature  at  the  crossing 
(interaction)  points.  Our  immediate  concern  then,  at  least  for  this  model,  is  to  de¬ 
termine  the  consequences  for  this  interpretation  of  our  transformation  to  a  solvable 
statistical  mechanics  model  and  the  conclusions  about  its  thermodynamics  that  we 
thereby  draw.  Our  results  should  be  compared  to  a  model  where  the  energy  of  a 
configuration  depends  on  the  cur\'ature  as  it  does  both  in  models  for  spacetime  and 
models  for  crystalline  materials.  Of  course,  this  also  raises  the  general  question  of 
what  energies/weights  might  be  appropriate  when  such  a  simple  lattice  structure  is 
not  available.  For  an  interacting  particle  model  of  a  generic  dynamical  system  the 
weakest  information  to  which  one  would  have  access  is  the  causal  ordering  of  the 
interactions;  our  work  to  this  point  suggests  addressing  the  question  of  whether  at 
least  some  such  models  could  be  mapped  into  a  statistical  mechanics  framework: 
though  one  could  not  expect  the  resulting  model  to  be  solvable,  some  understanding 
of  the  dynamical  system  might  be  gained  nevertheless. 


5.  References 

[1]  D.  A.  Meyer,  “State  models  for  link  invariants  from  the  classical  Lie  group^.'" 
in  A.  Kawauchi,  ed.,  Knots  90,  Proceedings  of  the  International  Conference  on 
Knot  Theory  and  Related  Topics,  Osaka,  Japan,  15-19  August  1990  (Berlin: 
W.  de  Gruyter)  559-592. 

[2]  S.  Takesue,  “Reversible  cellular  automata  and  statistical  mechanics. "  Pliys. 
Rev.  Lett.  59  (1987)  2499-2502. 

[3]  B.  Hasslacher  and  D.  A.  Meyer,  “Knot  invariants  and  cellular  automata.’' 


5 


Physica  D45  (1990)  328-344. 

[4]  B.  Hasslacher  and  D.  A.  Meyer,  ‘‘Knots,  criticality  and  thermodynamics  in 
discrete  dynamical  systems,"  submitted  to  Phys.  Lett.  A. 

[5]  G.  ‘t  Hooft,  “A  two-dimensional  model  with  discrete  general  coordinate-invari¬ 
ance,”  Duke  University  preprint  (19S9). 

[6]  B.  Hasslacher  and  D.  A.  Meyer.  “Lattice  gases  and  exactly  solvable  models," 
J.  Stat.  Phys.  68  (1992)  575-590. 

[7]  S.  Takesue,  “Ergodic  properties  and  thermodynamic  behavior  of  elementary 
reversible  cellular  automata.  I.  Basic  properties,”  J.  Stat.  Phys.  56  (19S9) 
371-402. 

[S]  S.  Wolfram,  “Statistical  mechanics  of  cellular  automata,”  Hev.  Mod.  Phys. 
55  (1983)  601-644. 

[9]  S.  Takesue,  “Fourier's  law  and  the  Green-Kubo  formula  in  a  cellular  automa¬ 
ton  model,”  Phys.  Rev.  Lett.  64  (1990)  252-255. 

[10]  R.  J.  Baxter,  Exactly  Solved  Models  in  Statistical  Mechanics  (New  "librk: 
Academic  Press  19S2). 

[11]  M.  Wadati,  T.  Deguchi  and  Y.  Akutsu,  “Exactly  solvable  models  and  knot 
theory,”  Phys.  Rep.  180  (19S9)  247-332. 

[12]  G.  E.  Andrews,  R.  J.  Baxter  and  P.  J.  Forrester,  “Eight-vertex  SOS  model 
and  generalized  Rogers-Ramanujan-type  identities,”  J.  Stat.  Phys.  35  (1934) 
193-266. 

[13]  L.  H.  Kauffman,  “Knots,  abstract  tensors  and  the  Yang-Baxter  equation,"  in 
L.  Lussana,  ed.,  Knots,  Topology  and  Quantum  Field  Theories,  Proceedings 
of  the  Johns  Hopkins  Workshop  on  Current  Problems  in  Particle  Theory  13 
(Singapore;  World  Scientific  1990)  179-334. 

[14]  M.  N.  Barber,  “Finite-size  scaling,”  in  C.  Domb  and  D.  L.  Lebowitz,  eds.. 
Phase  Transitions  and  Critical  Phenomena,  Vol.  8  (London:  Academic  Press 
1983)  Ch.  2; 

J.  L.  Cardy,  “Conformal  im-ariance,”  in  C.  Domb  and  D.  L.  Lebowitz.  eds., 
Phase  Transitions  and  Critical  Phenomena,  Vol.  11  (London;  Academic  Press 
1987)  Ch.  2. 


6 


